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ESTIMATION OF DROUGHT SEVERITY 

USING GIS 

II Hi 

1. Introduction 

Since the previous estimation of drought severity by Bangladesh Agricultural Research 
Council (BARC), which was published in 1990 in the title of "Drought in Bangladesh 
Agriculture and Irrigation Schedule for Major Crops" , there has been no attempts in 
updating the national drought maps. Since the initiation of the GIS Project, the climatic 
database such us rainfall, Evapo transpiration humidity, temperature, sunshine and cloud 
cover data collected from the respective organizations and thereby updated die previous 
climatic database, which was stored in the form of climatic resources inventory as part of 
the National Land Resources Inventory (LRI).. Under GIS project, attempt has been taken 
on updating of the previous drought maps incorporating latest climatic database and other 
dominant factors c.g. soil moisture-holding capacity, soil percolation rate. The mapping 
unit in the previous study was considered as soil association polygons accompaning the 
necessary attributes with overlaying climatic factors e.g P-Pre-kharif Transition period, 
K-Kharif and Rabi growing period, T-Cool temperature lines and c-cxtrcmc summer 
temperatures on to it. In the current GIS project, since all the soil and land resources 
database and also climatic database have successfully been transformed into GIS system 
and updated as of the present slate of data availability, this was the right time to update 
the drought maps analyzing latest climatic database and incorporating more layers under 
GIS platform. 

2. Objectives 

The objective of this research is to find out areas prone to drought in different seasons of 
the year. DruughL prone areas mainly located in the western part of the Bangladesh as 
estimated in the previous study. Drought is generally considered for the agricultural crops 
in three different seasons khraif-Il, Rabi and Pre-kharif. In the Kharif-11 (July-Octber), 
T.Aman is considered to be the major crops grown and is affected due to shortage of 
available soil moisture. Since T.Aman is grown in the wet season, soils having higher 
moisture holding capacity and low percolation rate with high to medium highland is 
mostly suitable for T.Aman and those arc normally not affected by drought. In the pre- 
k bar if season (March- June) Aus crops (Jute or Aus paddy) is affected. In Rabi season 
(Nov- Feb) Rabi crops (Whcat/Polato etc) arc normally affected due to shortage of 
available moisture in the soil. In the Rabi season, normally the whole period is non- 
rainfall situation, only the residual soil moisture supportAowingirops in this season. In 
the rabi season crops arc grown in such areas, where motfiure- holding capacity in the soil 
is higher and supplementary irrigation is ensured. 

Considering the above factors, GIS is selected to be used as a tool to analyze all the above 
mentioned factors and thereby delineate drought severity maps as follows: 

a) Kharif (T.Aman) drought prone area map; 

b) Rabi drought prone area map; and 

c) Prc- Kharif drought prone area map. 
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3. Scenarios Considered 


• Soil Characteristics that are mainly from the LRl database which is a survey results 
that reflects average condition in the soil. 

• Climatic Characteristics (LGP-Length of Growing period) that was calculated based 
on the time series database and calculated an arithmetic average that has 50% 
probability i.e average situation 

• Rainfall Uncertainty (rainfall availability at H0% probability), which represent 
ensured rainfall situation i c. this rainfall situation, is quite likely to occur almost 
every year. 

Based on the above scenarios discussed for different factors the drought severity is 
calculated and weight is given and accumulated in such a way that it represents an 
average situation that has the probability of occurrences is 50% i.e. it might occur once in 
two or two and half years. 


4. Factors considered 

In the current activity similar methodology followed on analyzing the factors as it was 
followed in the previous study in undertaking the tasks. The potentialities of using GIS in 
this task is that one can consider as many factors as one desires, analyze them anti this 
tools helps to delineate drought prone areas. 

In determining the potential factors experts from the project (B ARC’S core stalls and 
consultants) and national experts available in the country were consulted. As it is known 
that drought is a combined effect of climate and soil properties (specially factors those 
determines the moisture availability in soil). Therefore mainly the soil and climatic 
factors were considered in delineating drought prone areas. 


4 . 1 Kharif Drought area mapping 


Following factors have been considered in estimating the khar if 11 drought severity areas: 


Table I : Kharif ( I . A man) drought factors 


Soil Factors 

2 

Available moisture holding capacity 

-.idl- ’ — 

Soil percolation rale 

Climatic Factor 

3 

Kharif Length of growing period (days)-K Values 

Moisture Availability 
due to Rainfall 
uncertainly 

4 

Percentage of dry sub-humid (P<0.5PET) decades 

5 

Percentage dry decades (decades oT no rainfall) 
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4.2 Habi Drought area mapping 


billowing I aclo ns have been considered 
severity areas: 


in cslimaling the 


intensity of rabi drought 


I able 2 : Rabi drought factors 


Soil Factors 

l 

Available moisture holding capacity 

Climatic Factor 

2 

Rabi Length Length of growing period (days)-R values 
Percentage of dry sub humid (P<0.5PET) decades 

Moisture Availability 
due to Rainfall 
uncertainly 

¥ 

4 

Percentage dry decades (decades of no rainfall) 

External Water Supply 

5 

Ensured Irrigation 


3.3 Pre-Kharif Drought area mapping 

billowing lacio rs have been considered in estimating the pre-kharif drought severity 
areas: 


Table 3 : Prc-Kliarif drought factors 


Soil Factors 

1 

Available moisture holding capacity 

Climatic Factor 

1 2 

Prc- Kharil 'Transition period (days)-P values 


3 

Total Length of Growing Period (days)-LGP 

Moisture Availability 
due to Rainfall 

4 

Percentage of dry sub-humid (P<0.5PET) decades 

uncertainly 

5 

Percentage dry decades (decades of no rainfall) 


To consider the pre-kharif uncertainly of rainfall, the pre-kharif transition period (P) 
(numbers of uncertain days required to start full monsoon after the dry rabi period) is 
normally taken into account. Total LGP is considered' here, because the Total LOP 
determines the available atmospheric water that is potentially available for crops to be 
grown for the whole growing seasons. 

5 Estimate different parameters 


5. 1 Soil Parameter: Available Moisture Holding Capacity (AMHC) 


Using UR I summery program developed at the (JUS project, countrywide available 
moist lire holding capacity theme was produced. Five classes of moisture bolding capacity 
defined in the UR! as follows: 


Pa>:i' s of 17 


/ o 


M oisl u re A vailaliHity ra nge 
<HK) mm 
100-200 mm 
200-300 mm 
300-400 mm 
> 4(H) mm 

I he LR] database has one to many relationships from its mapping units with its soil 
attribute database due to fact that one soil association contains different percentage of 
occurrence of soil series. Therefore, for summarizing the attribute data to assgin a legend 
class to the map units, following category of dominance were defined such as a class 
occurring >80% is predominant class, a class occurring 60-80% is the mostly class, a 
class occurring >30% is with some class, a class occurring >15% is the mixed classes 
and a class occurring < 10 % is insignificance. Drought weights ( 1 - 5 ) signifying 
occurrence of different intensity of drought were also assigned, where I -very prone to 
drought i.e. low moisture holding capacity to 5- very low prone to drought i.c. high 
moisture holding capacity. 

The classes are obtained after summarizing is listed in Table 4. The corresponding weight 
to the classes that is assigned is also listed. An alternative way of assigning weight to 
different classes has been devised and applied also in this study. 

5.1.1 Revised mapping of the AMHC factor 

3 - An alternate way of assigning drought severity weight to different classes is devised 

^ based on the percentage of occurrence of different classes within each soil associations. 

The detailed table found after summarization (SUM2 fable) used as input and opened in 
excel, hive moisture classes were assigned to different drought severity weight as shown 
in Table 5. The procedure of mapping is explained in details in the six monthly reports of 
the National G1S consultant for Junc-Dcccmbcr ‘ 2000 . 

A cumulative weight was calculated under excel by multiplying percentage of occurrence 
of a class with its corresponding weight (Table 5). A part of the total worksheet is shown 
in Table 6 and brief description of the fields and calculation procedure is given at the 
bottom of the table. Weight assigned based on both the procedure is also compared in the 
Table. This procedure of applying weight is robust, logical and automated which can be 
used for other factors as well as for mapping of any theme of the LRI database, finally 
using the MHCWT New field the moisture holding capacity map is produced as shown 
in Figure 1. This theme was used as one of the input for drought severity 
cstimation/di ought mapping. 


Moist ure Cl asses 
1 

2 

3 

4 

5 
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1 able 4 : Moisture Classes for mapping based on Uic LR1 database 


CLASS 

WT 

CLASS 

WT 

Predominantly class 5 

5 

Mixed 5 2 

3 

Predominantly class 4 

4 

Mixed 4 3 2 

4 

IVcdominnntly class 3 

3 

Mixed 4 3 1 

4 

Predominantly class 2 

2 

' Mixed 4 3 

4 

Predominantly class 1 

1 

Mixed 4 2 3 

3 

Mostly class 5 

5 

Mixed 4 2 

3 

Mostly class 4 

4 

Mixed 3 5 4 

4 

Mostly class 3 

3 

Mixed 3 4 2 

3 

Mostly class 2 

2 

Mixed 3 4 

3 

Mostly class I 

1 

Mixed 3 2 5 

3 

Mostly 5 with some 3 

5 

Mixed 3 24 

3 

Mostly 5 with some 2 

S 

Mixed 3 2 I 

2 

Mostly 4 with some 3 

4 

Mixed 3 2 

3 

Mostly 4 with some 2 

4 

Mixed 3 1 2 

2 

Mostly 4 with some 1 

4 

Mixed 2 5 3 

3 

Mostly 3 with some 5 

3 

Mixed 2 4 3 

3 

Mostly 3 with some 4 

3 

Mixed 2 4 

3 

Mostly 3 will) some 2 

3 

Mixed 2 3 5 

3 

Mostly 3 with some l 

3 

Mixed 2 3 4 

3 

Mostly 2 with some 5 

2 

Mixed 2 3 J 

2 

Mostly 2 with some 4 

2 

Mixed 2 3 

2 

Mostly 2 with some 3 

2 

Mixed 2 1 3 

2 

Mostly 2 with some 1 

2 

Mixed 2 1 

2 

Mtwlly 1 with stxnc 4 

1 

Mixed 1 4 3 

2 

Mostly 1 with some 3 

1 

Mixed 1 3 2 

2 

Mostly 1 wifi some 2 

1 

Mixed 1 3 

2 

Mixed 5 3 2 

4 

Mixed I 2 3 

2 

Mixed 5 3 

4 

Mixed 1 2 

1 


4li 


Table 5 : Assigned weight to the moisture classes 


Moisture 

Classes 

Moisture 

Availability range 

Drought 

Weight 

Drought Severity 

1 

<1(M) mm 

I 

Very severe 

2 

100-200 mm 

2 

Severe 

3 

2(X)-3(X) mm 

3 

Moderate 

4 

300-400 mm 

4 

Low 

5 

> 400 mm 

5 

Very Low 
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5.2 Soil Parameters: Soil Percolation rate 


¥ 


Sml [u uoliiliou i ale is anolhei important factor to estimate drought because normally drought 
m the T.Aman crop period (KlmriO occurs in the high permeable soil (sand, sandy 
loam/sand y silt etc.) rather than low permeable soils (clay, silty clay etc), where other 
conditions remains the same. In the LRI, information on soil percolation is given as the 
heading of soil permeability; three classes of permeability arc defined in the LRI c.g. slow 
(<i: cm /day), moderate (12-305 cm/day) and rapid (305 cm/day). It is understood that these 
classes are very lumped and too much generalized to use this information for drought severity 
estimation therefore through consultation with soil scientist and based on different 
literature, a more detailed percolation rate theme was worked out based on the topsoil texture. 
Live different percolation rate were identified based on the 9 different top soil texture 
property those were considered to be used in this study. The corresponding weights to soil 
permeability class were also assigned a weight factor for drought prone areas due to the soil 
percolation rate. The relationship between soil textures and soil percolation arte, drought 
severity weights arc shown in Table 7. 


fable 7 : Soil properties in relation with the soil texture class 



Coll 

Col- 2 

Cob 3 

Col-4 

Col-5 

Col-6 

' SL 

Texture 

Field 

Willing 

Available 

Percolation 

Drought ratings** 

No 

( ’lass 

Capacity* 

Point* 

Capacity 

* 

Rate** 

(Drought prone 
areas) 

7 

Sand 

0.12 

0.04 

0.08 

Very high 

1 

0 

A* 

Loamy Sand 

0.14 

0.06 

0.08 

High 

2 

j ” 

Sandy Loam 

0.23 

0.10 

0.13 

High 

2 

4 

Loam 

0.26 

0.12 

0.14 

Medium 

3 

5 

Silt Loam 

0 30 

0.15 

0.15 

Medium 

3 

6 

Silt 

0.32 

0.15 

0. 1 7 

Low 

4 

7 

Silty Clay 
Loam 

0.34 

0.19 

0.15 

Low 

4 

X 

Silty Clay 

0.36 ~ 

0.21 

0.15 

Very Low 

5 

0 

Clay 

0,36 

0.21 

0.15 

Very Low 

5 


*Vnhtr.s obtained from ASCK /9WJ, Table 2 6, p.2, ** Values assigned in consultation with the soil scientist 


In the LRI there are 20 different texture classes are defined and 9 classes o I textures as shown 
in Table 7 arc related to 20 textural classes as shown in Table 8. 


Table 8 : LRI Soil texture class and corresponding texture class of Table 7 


AEZ 

Class 

Description 

Texture 
class in 
Table 7 

AEZ 

Class 

Description 

Texture 
class in 

Table 7 

1 

n 

3 
1 4 
s 

Sand 

1 

11 

Sill loam 

5 

Loamy sand 
Loamy fine sand 
Sandy loam 
Line sandy loain 

2 

2 

3 

3 

12 

13 

14 

15 

Gravelly clay loam 


Silty clay loam 
Clay loam 
Silty clay 

7 

7 

X 
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AEZ 

Class 

Description 

’Texture 
class in 
Table 7 

AEZ 

Class 

Description 

Texture 
class in 

Table 7 

6 

Very line sandy 
loam 

3 

li 

Clay 

9 

7 

Gravelly sandy 

clay loam 

3 

17 

Mucky clay 

9 

8 

Sandy clay loam 

h 3 

18 

Muck 

9 

9 

Loam 

4 

19 

Peaty Muck 

9 

10 

Sill 

6 

20 

Peal 

9 


Following steps were followed to create percolation rate map from the topsoil texture data of 

LRI. 

• Use LRI Summary program, select TX-Topsoil texture as a classification field, and 
choose AEZ-soil association, as the mapping polygon and maps and tables are die output 
options. 

• Use assign reclassification (regroup) option (c.g. Yes) and assign values as per Table X 

• Use default cutoff percentages 

• In the editing class dialog box, choose “Continue" and at the end two tables 
NTXJVEZCODESUM2 and NT.X_AEZCodcSUM3 tables will be displayed and a view 
named NTX AHZCode View will fie displayed, 

• Use NTX AEZ('odeSum2 Table and use Ihc alternative procedure (explained in section 
*1.1 and Table 6). 

• Reclassify the 9 texture classes with die 5- percolation rale classes as given in table 7, 
col 5. 

• Finally the percolation map is obtained by assigning the percolation rate class applying 
the same procedure ot weight application to different percentage of occurrence and 
thereby finally assigning the legend to the map unit as explained in Table 6. 


The percolation rate map is shown in Figure 2. This percolation rale theme was used as one 
of the theme for analysis in delineating the drought prone areas. 


For further documentation on the LRI Summary program and steps, please see six monthly 
teporl of the System Analysl/Lrogrammer of July-Dcccmber 20(H). 

In summary weight grids of the soil themes arc calculated based on the weight ol the soil 
available moisture holding capacity and soil percolation rale. J heme as shown in Table 9. 
These weight grids are used later on for computing the drought severity areas. 


Table 9: Assignment of weight to different soil factors 


Factor 

Weight 

Soil Moisture 


Soil Percolation talc 


- 

Values of different fac 

tors 




1 

2 

3 

4 

5 


< 1 f X) mm 

1(H) 200 
mm 

200-300 

mm 

300-400 

mm 

>4(H) mm 

e 

Very high 

111 ® . 

Moderate 

Low 

Very Low 


farter 
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5.3 Climatic Parameters 


5,3.1 Kharif Length of Growing Period 

Analyzing climatic data of I960 1990, kharif length of growing period was calculated lor all 
BWDB and BMD .stations (291 stations) using API 4.0 software (developed hy FAO lor 
calculation of Net Biomass, EvaptMranspiralion, Crop Water requirement etc.). Climatic data 
in APT4.0 format was provided and using the CDA module of APT4.0 kharif LGP was 
calculated. Average kharif LGP is shown in Appendix- A. The average kharif LGP Map is 
shown in Figure 3. In Figure 3, there are 12 different classes of K-LGP zones (each having 15 
days interval) are round from the data. Considering the fact that as higher the LGP lower the 
drought probability, a drought severity weight to the K-LGP values is assigned as shown in 
Table 10. 


Table 10: Weight Assigned to Kharif LGP surface 


SL 

No 

Kharif LGP 
(Days) 

Weight 

SI No 

Kharif LGP 
(Days) 

Weight 

\ 

185-200 

1 

7 

275-290 

4 

2 

2QO-2 1 5 

l 

8 

290-305 — 

4 

3 

215 230 

2 

9 

305-320 

5 

4 

230-245 

2 

10 

320-335 

5 

| 

245-260 

3 

1 1 

335-350 

5 

6 

260-275 

3 

12 

350-365 

5 


I -Very prone to drought, 2 Moderate drought, 3-Lcss Drought, 4-Slight 
5-No drought. 


5.3.? Total Length of Growing Period (Rabi + Kharif) 


As mentioned earlier that this factors was considered for estimating pre-kharif drought 
severity estimation. Based on historic climatic data analysis (1960-90) average total length of 
growing period (days) calculated using APT4.0 and the corresponding LGP surface is shown 
in Figure 4. This factor signifies the overall moisture situation that could be used in prc-khnrf 
period. 

5.3.3 Rabi LGP and corresponding weight to drought severity 

In Figure 5, there arc 5 different classes of Rabi LGP zones (each having 10 days intervals) 
are found from the data. Considering the fact that as higher the Rabi LGP lowers the drought 
probability in the rabi season, therefore a drought severity weight to the Rabi LGP- values is 
assigned as shown in Table I I. 

Table 1 1: Weight Assigned to Rabi LGP surface 


SI. 

No 

Rabi LGP 
(Days) 

Weight 

Si No 

Rabi LGP 
(Days) 

Weight 

1 

80-90 

T 

5 

120-130 

5 

2 

90-100 

2 

6 

130440 

5 

3 

100 1 10 

3 

7 

140-130 

5 

4 

1 It) 120 

t 




L_ — - — — - 


I .Very prone to (irnugi)i. 2 -Moderate drought, 3 Less Drought, 4 Slight, 5-No drought 
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5.3.4 Dry Sub-humid decades and Moist Sub-humid decades 


The dry sub-humid decodes mid moist sub- humid decodes is different cropping seasons 
factors determine the stress situation that could be prevailing in the corresponding session 
concerned. Decadal calculation is the modest interval that has been considered for the 
analysis in the current study. An atmospheric moisture balance is assumed as shown in 
Figure-6 to calculate the above parameters. In Figure 6, P am! PET vs. time arc plotted. The 
days having P < 0.5 PET is considered as dry Sub- humid period. The days having P > 0.5 PEI' 
and P <= PET is considered as moist sub- humid period and the no of decades having no 
rainfall is called dry decades. The number of days having P > PET is called humid period. 
Average monthly PET (calculated using Penman- Montielh formulas) has been calculated by 
APT4.0 using daily time series climatic data of 1961-1990. Decadal PET is calculated the 
division by 3 of the monthly PET value. The P value has been considered here as the 

dependable rainfall that is the rainfall at 80% probability. 

i 



Figure -6: Conceptual diagram of atmospheric moisture balance 

5.3.4. 1 Decadal Rainfall and Decadal PET Calculation 

From the decadal rainfall data of different years of record, rainfall at any frequency of 
occurrence for as many as stations arc desired is expected has been calculated using a 
FORTRAN program (Appendix A). Following arc the necessary steps to calculate rainfall 
frequency file and a decadal ET value from the decadal rainfall and monthly PET data. 
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Prepare inputs for the program: Average monthly PET file, a listing of rainfall stations Hies 
ami decadal rainfall tiles (as per AIM 4.0 decade file format) are the inputs. 

Run the rfpdce.exe. give the input files, first ET file and secondly the file having the rainfall 
file listing 


Note: Input files should he located at same directory location, from where the program is 
being executed 


(live Output decadal PET file and output Rainfall flic name and then it asks the percentage of 
probability is expected Finally it gives expected results. 

Two output files are calculated one is the decadal ET and another is decadal Rainfall at the 
required probability. One extra column is added in the rainfall file, which is number of year 
records. If number of years record is >= 30 years, the probability values are significant. 

An example of decadal PET and rainfall at desired probability calculation: 

For example to calculate decadal rainfall with 75% probability and decadal PET of the 
BWDB stations 


Following are the inputs: 

it Hie BWDB rainfall data II Ic is stored in NT2\\C:\Hasan\Clim-An:iIysis\HWDB 
all the dee files are the historical rainfall data in APT4.0 Decadal formal 

ii) Rnlist is the Rainfall file list 

iii) BWDB-ET.pni is the monthly ET file 

iv) loiter output decadal ET "wdbcl.dee" and output rainfall file "wdbrn75.dee“ and 
type 75, when it asks the probability of occurrences. 

Following are the outputs: 

I wo files ax inputted catlicr "wdbct.dec” and "\vdbrn75.dec" arc eretaed during execution of 
the program 

I he " w dlict. dec' * mil put file is as shown below : 

l . n if, q 35.9 3 2 6 32.6 Up l o 36 columns 

^ 

ami the M wdbrn75.dcc" output file is as shown below: 

j 13 n ,o .0 .0 .0 .0 Upto 36 co him 

it; n n .0 .0 .0 .0 Up to 36 colum 

* 

*1 llltfl »•* * 

1. 
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Similar calculated is to tic followed for the BMD station files are stored in 
NT2\\C:\I Iasan\C'lim-Analysis\BMD directory 


5.3.4. 2 ( nictitation of Decadal Water balance 


I tic hies created above were used to calculate lire decadal atmospheric water balance data A 
l-ORTRAN program was written to calculate P-0.5 PET, where P = Dependable rainfall ami 
IT. I is the decal PE I'. Listing of the program is given in Appendix- B. It takes the output files 
generated as discussed in section 4.2.2. 1 as inputs and creates the a output Tile containing four 
rows of information for each stations, First row gives the rainfall (at the requested 
probability), second row gives the decadal PET values, third row gives the 0.5PET values and 
the fourth row gives the P-0. 5 PET values. 


as shown in below : 


Output listing is 

11 .0 . 0 

12 25.9 25.9 

13 12.9 12.9 

14 12.9 12.9 


.0 .0 .0 

25.9 32.6 32-6 

12.9 16.3 16.3 

-12.9 16.3 -16.3 


. 0 

.0 

.0 

32.6 

46.4 

46.4 

16.3 

23 .2 

23 .2 

-16.3 

-23 .2 

-23.2 


10.0 . .p Values 
4 6.4 . . PET va 1 ups 

0.5PET values 

- . P- 0 . 5 PET values 


This decadal water balance information can be mapped using GIS and also can be presented 
in graphs. In graphs time scries (year round) water balance is shown (Figure 7). In map form 
spatial variation of water balance condition at a particular time (e g. at any decade) can l>e 
presented.. Both arc useful in assessing the atmospheric water balance condition at a 
particular location for the specified period. 



figure 7: An example of water balance plot of I he Barisal Station at 50% rainfall probability 
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5. 3.4.3 Calculation of percentage dry decades and percentage dry sub-humid decades 
Using the output (lie of the decadal water balance (section 4.2.2. 2). frequency of (number) 
dry decades, dry sub-humid decades, moist sub-humid decades and humid decades arc 
calculated. A program called droughlpar.exe is written for this purpose. The source code is 
given in Appendix-C. The input, output files and the cxc file resides at NT2\\C:\Hasan\Clim- 
Ana!ysis\Droughtpar\debug directory. During execution of the program, the user is asked to 
provide start decades and end decades for calculation in addition to the output file containing 
i he number of decades having dry sub-humid, dry decades and most sub humid and humid 
decades. Moist sub- hum id and humid decades arc not used in I he drought severity 
calculation. 


An Example of calculation procedure: 


Say for example drought parameters for the BWDB stations for the pre-kharif transition is to 
be calculated. At first it is to In' ensured that the program file and the input Tile reside at same 
directory location. In this case the files resides at NT2\\C:\Hasan\Clim- 
Analysis\Droughlpar\dcbug directory. Now run droughtpar.exe from the DOS prompt or 
from the file explorer window. During execution the program asks for the name of the input 
file and the output file. Enter input Tile as BWDBWB80 (where decadal water balance 
information are found, as discussed in 4.2. 2.2) and output file as PK0715.txt, enter 7 as start 
decades (start of March) and 15 as (he end decades (End of May), ‘flic output file consists of 
live columns as shown below: 


10 18 5 

2 0 3 6 5 

3 0 0 9 6 

4 0 18 4 

5 12 6 4 

6 0 2 7 5 

7 0 0 9 5 

9 0 3 6 4 

10 0 2 7 4 


0*1 1 : .Station axle ; 0*1 2 Number ol Humid decades. 0*1 3 : Number of Moist Sub-Humid, 0*1-1 :Numl*er 
i*r Dry sub humid and 0*1 5 Number ol dry decades. 

The first column (station ID) is used to link with the corresponding shape file or coverage lor 
surface creation of percentage dry decades and percentage dry sub-humid decades. 
Percentages are calculated after opening this file in EXCEL program. 


It is to l>e noted that the sum of the values of humid, moist sub-humid and dry sub humid is 
equal to (he total period, because the dry sub-humid decades also includes dry decades. I be 
fifth column (dry period) is the total no rainfall decades within the analysis period 


Using the above data file and the station shape file, percentage dry decade and percentage dry 
sub humid decade surface for both pre-kharif and kharil seasons were prepared as shown in 
figure 8. 9, 10 and I 1 Percentage of dry decades and percentage of dry .sub humid decades 
:dso calculated lor the Kabi period but since most of the time of Ibis season (i.c, October 
Eehruar v) is nmmally no rainfall, this factors carries less importance than the raid season 
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trrigiti ion. The rahi season irrigation percentage based on the NM1C (National minor 
irrigation census) data of 1996/97 season is shown in the map (Figure 12). This was also 
considered during the Rahi season drought estimation. 


Percentage of dry sub- humid decades and dry decades was considered to quantify drought 
severity for the season concerned. Percentage cut off points of the dry sub-humid and dry 
decades classes for quantifying the drought severity is presented in Table 12. 


Table 12 : Percentage cut off for the drought classification 


Drought 

weight 

Drought prone areas 

Percentage cutoff of the dry factors 

Dry sub-humid 
period 

Dry decades 

i 

Very severe drought prone 
areas 

>80 

>80 

2 

Severe drought prone areas 

65-80 

65-80 ~~ 

3 

Moderate drought prone 
areas 

45-65 

45-65 

jT 

Slight drought aicas 

30-45 

30-45 

5 

No drought prone areas 

<30 

<30 


Using the above cutoff points six weight surfaces (two for each season, one for dry sub 
humid and one for dry decades), showing drought prone /.one due to the above factors for 
each seasons was calculated. 

In a concise form the weight value for all the climatic parameters calculated using the above 
methodology is given in Table 13. 

Table 13: Assignment of weight to individual climatic factors 


Factor 

Values of different factors 

Weight 

1 

2 

3 

4 

5 

Pre-Kharif Transition 
period 

>60 days 

50-60 

days 

40-50 

days 

30-40 

days 

< 30 days 

Kharif LGP 

”<200 

days 

200-220 

days 

220-240 

days 

240-260 

days 

> 260 
days 

Rahi LGP 

<90 days 

90- 100 
days 

100 1 10 

days - 

I 10-120 
days 

> 120 
days_ 

Total LGP 

240-255 

days 

255-270 

days 

270-285 

days 

285-300 

days 

> 300 
days 

Pre-Kharif Dry Decades 

>71 % 

63-71 % 

46-63% - 

29-46% 

<40% 

Pre- Kharil Dry Sub- 
humid 

>90% 

80 90% 

60-80% 

40 60% 

<40% 

Kharif Dry Decades 
Kharil Dry Sub-humid 


15-18% 

12-15% 

9-12% 

<9% 

22-30% 

19-22% 

15-19% 

12-15% 

. 4 . . .i — r. 

<12% 

. ( 1* n fMfc J If* t f"i 


f;i« I t»r 
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6. Step by step procedure for preparation of drought maps. 

6. 1 Preparation of Pre-Kharif Drought Map 

The time period for pre-khar if drought is 2 ml decade of March through 1st decade of June 
(Week 8 to Week 16) period. 

• Calculate percentage dry sub-humid and percentage dry decades surface using the 
procedure mentioned in section 4 and classify the surface grid according to weight 
signifying severity of drought as mentioned in Table 13. 

• Calculate Soil-association wise moisture holding capacity surface and classify 
according to drought severity weights (Table . 

• Calculate Prc-kharif growing period surface and assign weight to different classes 
according to Table 13. 

• Calculate overall prc-kharif drought index, by combining the above mentioned four 
weight grid surfaces using the following formula 

(MllCWt * 40) 4 {KHLGPWi * 25) 4 {DSUWt * 20) 4 (DDWt + 1 5) 
PKHDrlndex - - — 

Minimum value for the prc-kharif overall index would be 10 considering very severe drought 
situation prevailed from all the factors concerned (weight- 1) and maximum value would lx' 
50 considering very good moisture situation for the crop (wcighl-5). 

An overall prc-kharif drought index is obtained starting values from 16 to 50, where a 
severity is assigned in terms of 5 different classes of the above indices. 1 6-20- Very severe, 
21-26 Very severe class, 26-31 Moderate Class, 32-40 Slight class and 41-50 is no drought 

class. 

The final Pre-Kharif drought map found from the above calculation is shown in Figure 13. 


6.2 Preparation of Kharif Drought Map (T.Aman) 


The time period for Kharif drought is 2 nd decade of June- end of October for the T.Aman 
Crop (decade 17 to decade 30). Following different surfaces arc calculated: 


Calculate Percent dry sub-humid and Percent dry decades surface and classify the 
values according to weight signifying severity of drought as mentioned in Table 13. 

Ilif '7’ : F 

Calculate Kharif length of growing period surface and classify according to severity 

of drought (Table 13) , r 

Calculate Suil-assoeialion wise moisture holding capacity surface and classify 

according to severity (Table 9) 

Combine the above mentioned four surfaces to create a cumulative drought index 
using following formula. 
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I'KHDrlmle t = (A, /,( lVl * 40 > + ( KIll.GI’Wl * 25) + (l)SIIWl * 20) + U)l)Wt * 1 5) 

10 

Reclassify the cumulative grid according to the level of severity (10 16, Very severe, 17-21 
Severe, 22-27 Moderate, 27-34 Slight and >34 is nil or no drought.) 

In this map some not T.Aman areas were excluded from the mapping considering the fact that 
this map is prepared to show drought severity to T.Aman crop only. Not T.Aman areas arc 
those areas having landtype Lowland to very low land. Hilly regions arc also excluded from 
the mapping. Therefore a non-T.Aman area grid is obtained first and after doing some grid 
operation ( mainly the setmill and con request) the final T.Aman drought map is obtained as 
shown in figure 14. 


6.3 Preparation of Ftabi Drought Map 

The time period lor Kharif drought is 1 sl decade November- I s ' decade of March for the 
major rabi crops (decade 31 to decade 7). following different surfaces arc calculated: 

• Calculate Percent dry sub- humid and Percent dry decades surface and classify the 
values according to weight signifying severity of drought as mentioned in fable 13.. 

• Calculate Rabi length of growing period surface and classify according to severity of 
drought (Table 13) 

• Calculate Soil-association wise moisture holding capacity surface and classify 
according to severity (Table 9) 

• Combine the above mentioned four surfaces to create a cumulative drought index 
using following formula. 

(MHCWt * 40) + {KHLCrWt * 25) + (/ )SHWt * 20) + {l)DWt * 1 5) 
PKHDrlnd rX ~ ” 


Reclassify the cumulative grid according to the level of severity (10-16, Very severe, 17-21 
Severe, 22-27 Moderate, 27-34 Slight and >34 is nil or no drought. 


finally the drought prone area map for the Rabi season is obtained as shown in figure 15. 
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Appendix A 

Program listing for Decadal Rainfall and Decadal ET creation 

PROGRAM RNPROB.FOR 
C 

Character* il ETC i ie, RN1 ist , Dummy, ETout , RNout 

Character *9 Rnfile 

INTEGER Dec (12, 3 ) , NSTP (3 00) 

DIMENSION D PET (300,36) , PET (300, 12} , DRAIN ( 40 , 3 6 ) , PRAIN ( 4 0 , 3 6 ) , 

+ DVAL (36) 

write(*,*} 'Enter the ET-File 
Read ( * , M ETEile 

write (*,*)' Enter Rainfall file List' 

Read ( * , * ) RNlis t 

writet *, *)ETf ile, RNlist 

open (3 , FILE-ETf i le, STATUS^ 'OLD' ) 

open (4 , FILE=RN1 ist, STATUS='OLD' ) 

write{*,*) 'Enter output decadal ET-File 

Read { * , * ) ETout 

write (*,*)' Enter output decadal Rainfall file 
Read ( * , *) RNout 

wr i te {*,*)' Enter probability (%) of decadal Rainfall 
Read ( * , * ) RNProb 

Open (6, FILE = RNout, STATUS= ' NEW' ) 
open (7 , FILE-ETout , STATUS^ 'NEW ' ) 

C Calculate Decadal PET from the Monthly PET and then write it 

read (3 , * ) Dummy 
1=1 

90 Read (3, 120 , END-130 ) NSTP ( I ) , ( PET ( I , J > ,J=1,12) 
write{*,120) NSTP { I ) , (PET(I,J> , J=l, 12) 

1 = 1 + 1 

120 format (IS, 12F6. 1) 

Goto 90 

C Create Decadal file j, 

130 NST a 1-1 

DO 160 1=1, NST 
NR=0 

DO 160 J = l, 12 
DECPET=PET(T , J) /3 . 0 
Write ( * , * ) DEC PET 
C NC = 0 

Clll NR=NR+ 1 

C NC=NC+ 1 

CC DPET { I , NR) = DEC PEI’ 

C IF (NC . LT . 3 ) GOto 111 

NC=J * 3 
NNC=NC-2 
DO 140 N=NNC,NC 
DPET (I ,N) = DEC PET 
140 continue 

c write ( *, *)NSTP( I) , (DPET(I.M) ,M=1 ,36) 

160 continue 

C Writing the decadal PET data 

po 190 1=1, NST 


Wr i 1 n ( 7 , 180) NSTP Cl), (DPET (1,0) , J=1 ,36) 
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180 format (15, 36F6.1) 
190 continue 
Close ( 7 ) 


C ******** End 0 f the Decadal ET Calculation and writing 


C ************ start of Rainfall Data reading, sorting and calculating prob 
C NPST=0 

5 read (4 , 10 , END=1 000) Rnf i le 
10 Format |A9) 

C NPST=NPST+1 

open ( 5 , F I LE - R n f i 1 e , STATUS = ' OLD ' ) 

C Initialize the Sums 

DO 50 I “1 , 36 
DVAL ( I ) - 0 
SO CONTINUE 


NYR -0 

100 NYR- NYR + 1 

DO 200 1 = 1 , 12 

READ (5,20, END= 999) NST, (Dec(I,J) ,J=1,3) 
20 Format (4X, 15, 8X, 317) 


200 CONTINUE 
NR -0 

IX) 500 1-1,12 
Do 500 J=1 , 3 
NR -NR ♦ 1 

DRAIN (NYR, NR) -Dec (I , J) ! DRAIN ( 3 00 , 3 6 ) -> Decadal Rain 

S00 CONTINUE 


GOTO 100 

999 NYR - NYR -1 

C * * * End of reading one station data where NYR - No of years, NR- 
Decadess 

C Sorting tlie Data for Probability Calculation PRAIN ( NYR , NR ) - 

> DRAIN (NR ) 

DO 350 NR = 1,36 

DO 330 I = 1 , NYR 
Max -DRAIN ( I , NR) 

Maxpos = I 


300 

320 

C 

i ro 


K= I *1 

DO 320 J-K , NYR 
IF (DRAIN (J, NR) . GE .Max) THEN 
Max = DRAIN fJ, NR) 

Maxpos - J 


ELSE 

GOTO 300 
END IF 
CONTINUE 
CONTI NUE 

PRAIN (T ,NR) -Max 
SWAP the data with the 
(■ omp - DR A T N ( f . NR ) 

DR A I N ( Maxp* > s , NR ) “ 1 emp 
CONI’ I NUE 


!PRAIN-> is t)ie max of the serieses 
Initial position & the Max position 
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c 


Calculate rank, with the given Probability value and extract the data 
RANK -RN I 'ROB* (NYR 1 ) /I 00 
3 RANK - I NT (RANK) 

DVAI. (NR) = PRA IN ( 1 RANK , NR ) 

3 SO CONTINUE 

WRTTE(6, 610) NST, NYR, (OVAL ( I ) ,1-1,36) 

610 Format (15, IX, T3, 36F7.1) 

close ( S) 

GOTO 5 

1000 Close (4) 

Cl OSP { 6 ) 

Stop 

End 
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Program listing for Decadal Water balance calculation 


Appendix B 


Program PETDATA . FOR 

Character * 1 1 ETBMD , RNBMD , ETBWDB , EBWDB , RNOUT 
Character * 9 BMDDEC , BWDBDEC 
INTEGER NSTP(300) ,NSTR(300) 

REAL PET (300,36}, RAIN ( 3 DO , 36 ) , RNMPET ( 300 , 3 6 ) , UPET (3 00 , 3 6 ) 
write!*,*) 'Enter the decadal ET-File 
Read ( * , * ) ETBMD 

wr ite (* r *)' Enter decadal Rainfall file 
Read ( * , * ) RNBMD 

write (*,*)' Enter output file file 
Read ( * , * ) RNOUT 

c write! *,*) 'Enter the BWDB decadal ET-File : ' 

C Read < * , * ) ETBWDB \ 

c write!*, *) 'Enter BWDB Rainfall file :' 

c Read ( * , * ) RBWDB 

open {3, F1LE= ETBMD, STATUS = ' OLD' ) 
open (4 , FILE= RNBMD, STATUS='OLD' ) 
c open ( 5 , F I LE = ETBWDB , STATUS = ' OLD ' ) 

c open ( 6 , FILE=RBWDB, STATUS^ ' OLD' ) 

open (7 , FILE=RNOUT, STATU S= ' NEW ' ) 

C Read the BMD file first in the memory variables 

NSTPET-0 
DO 20 1=1 ,300 

read (3 , 10, END=30)NSTP(I) , (PET(I,J) ,J=1 ,36) 

10 Format ( 15, 36F6 . 1 ) 

20 Continue 

30 NSTPET=I-1 

NSTRN=0 
DO 50 1=1,300 

read (4 , 40, END=60) NSTR ( I) , ( RAIN ( I , J ) , J=1 , 3 6 ) 

40 Format ( 15 , 4X, 36F7 . 1 ) 

50 Continue 

60 NSTRN= I - 1 

Print*, 'No of BMD ET Stations NSTPET 
Print*, 'No of BMD RN Stations ,NSTR 
DO 100 1=1, NSTPET 
DO 110 J- 1 , NSTRN 

IF(NSTP(I) .NE.NSTR(J) )goto 110 
write (7, 120) NSTR(J), (RAIN { J , K) , K=1 , 36) 
wr i te (7 , 130) NSTP(I), ( PET ( J , K) , K=1 , 3 6) 
do 90 K~1 , 3 6 

RNMPET! I, K) =RA I N ( J , K ) -0 . 5 *PET { I , K) 

HPET ( 1 , K ) =0.5* PET ( I , K ) 

90 continue 

write (7 , 150) NSTP ( I ) , (HPET { I , K) , Kd , 36) 
write (7, 140) NSTP ( I ) , (RNMPET ( I , K) , K=1 , 3 6) 


110 

100 

120 

130 

150 

140 


goto 100 
continue 
continue 
format ( 15 , ' 
format (IS,' 
format ( 15 , ' 
f ormat (IS, 
STOP 
END 


1 

2 

3 


, 3 6F6 . 1 ) 
, 3 6F6 - 1 ) 
, 36F6.1) 


4 ' . 3 6F6 . 1 ) 
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Program listing for Drought parameter calculation 


Appendix C 


c Program DROUGHTPAR . FOR 

c Takes input from the WHBO.Out file, there are four rows data exists 

for all the stations 

Character* 10 InFi le, OutFile, Dummy 
REAL P(36), PET (36) ,1 1 PET (36) 

Write(* , * ) 'Enter the Input File Name:' 

Read ( *,11) InFile 
11 FORMAT (A10) 

Write (*,*)' Enter Analysis Output File Name 
Read {* , 11) OutFile 
Open (3 , FILE- InFi le, STATUS^' OLD' ) 

Open ( 4 , FILE-OutFi le, STATUS^ ' New ' ) 

Open ( 5 , F I LE= ' PRN ' ) 

Write (5, *)' Station Humid MSH DSH Dry Days' 

Write (*,*) 'Analysis Basis are following (No of Days):' 

Write!*,*) 'Humid -> P>PET' 

Write (*,*)' Moist Sub-Humid -> P>0.5PET' 

Write (*,*)' Dry Sub-Humid -> P<=0.5PET' 

Write(*, *) 'Dry Days -> P=0' 

Write t *,*) 'Mon : Jan Feb Mar Apr May Jim Ju) Aug Sep Oct Nov Dec' 
Write!*, *) ' Dec : 1 3 4 6 7 9 10- 13- 16- 19- 22- 25- 28- 31- 34-' 

Wr i te (*,*)' Pre Kharif & Rabi Drought -> March-May (Normal)' 

Write (*,*)' Khar if Drought occurs during ->Sept-Nov (Normal)' 

c Enter the period (Start Decades and End Decades) for which analysis is 
to be done . 

Write!*, *) 'Enter Start Decades Cor Analysis 
read ( * , *)NSTD 

Write!*,*) 'Enter End Decades for Analysis 
read ( * , * ) NEND 

C Write!*,*) 'Output file contains STNNo, Humid, 

Mo i s t Subllumid , Dr ySubHumd , Dry Decades ' 


100 

10 

15 


21 

20 


Write!*,*) 'Reading input file,.' 

Read (3, 1 0 , END-^999 ) NSTP, (P ( I ) , 1=1 , 3 6) 

Read (3 , 1 0 , END-999 ) NSTP, {PET{I) , 1=1,36) 
READ (3, 10, END-999) NSTP, (HPET(I) ,1=1,36) 
Format ( 1 5 , 2 X , 3 6F6 . 1 ) 

READ (3,15, END =9 9 9 ) Dummy 
Forma t ( A10 ) 


NHUM=0 
MSH = 0 
NDSH=0 
NDD=0 

DO 20 I = NSTP, NEND 
IF(P(I) .GT.PET(I) ) GOTO 21 
I F ( P ( I ) . EQ * 0 . 0 ) NDD=NDD+ 1 
IF(P(I) .GT.HPET! 1) )MSH-MSH< 1 
T F ( P ( T) . LE . HPET ( I ) ) NDSH=NDSH t 1 
GOTO 20 



mtUM=MHUMt 1 
CONTINUE 

WR 1 TE ( * , * ) NSTP , NHUM , MSH . NPSH , NDD 
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40 

999 


WR ITE (5 , * ) NSTP , NHUM , MSH , NDSH , NDD 
WRITE ( 4 , 4 0 ) NSTP , NHUM , MSII , NPSll , NDD 
Format (15,414) 

GOTO 100 

WR ITE ( * , * } '*** END OF PROCESSING 

STOP 
END 






* * * * 


